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Abstract. XTE J1810— 197 and IE 1547.0—5408 are two transient AXPs exhibiting radio emission with unusual 
' properties. In addition, their spin down rates during outburst show opposite trends, which so far has no ex- 

, planation. Here, we extend our quark-nova model for AXPs to include transient AXPs, in which the outbursts 

^\ • are caused by transient accretion events from a Keplerian (iron-rich) degenerate ring. For a ring with inner and 

outer radii of 23.5 km and 26.5 km, respectively, our model gives a good fit to the observed X-ray outburst from 
XTE J1810— 197 and the behavior of temperature, luminosity, and area of the two X-ray blackbodies with time. 
The two blackbodies in our model are related to a heat front (i.e. Bohm diffusion front) propagating along the 
ring's surface and an accretion hot spot on the quark star surface. Radio pulsations in our model are caused by 
I ■ dissipation at the light cylinder of magnetic bubbles, produced near the ring during the X-ray outburst. The delay 

O ' between X-ray peak emission and radio emission in our model is related to the propagation time of these bubbles 

to the light cylinder and scale with the period as tprop. oc P^~'^ j P^l"^ where a defines the radial dependence 
of matter density in the magnetosphere (oc r""); for an equatorial wind, a = 1, we predict a ~ 1 year and ~ 1 
month delay for XTE J1810— 197 and IE 1547.0—5408, respectively. The observed flat spectrum, erratic pulse 
profile, and the pulse duration are all explained in our model as a result of X-point reconnection events induced 
by the dissipation of the bubbles at the light cylinder. The spin down rate of the central quark star can either 
increase or decrease depending on how the radial drift velocity of the magnetic islands changes with distance from 
the central star. We suggest an evolutionary connection between transient AXPs and typical AXPs in our model. 
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QQ 1. Introduction Ouyed et al. 2002) building on three previous papers where 

, we explore its application to Soft Gamma-ray Repeaters 

^ , Anomalous X-ray pulsars (AXPs) are magnetars with ro- ^^^^^ (q^^^^^ Leahy, & Niebergal 2007a; OLNI), to 

tation period of 2-12 seconds and inferred surface mag- ^^Ps (Ouyed, Leahy, & Niebergal 2007b; OLNII), and to 

^ . netic field strength B ~ 10" G (e.g. Woods & Rotating Radio Transients (RRATs) (Ouyed et al. 2009; 

^ ; Thompson 2006; Kaspi 2007). In this work we focus on OLNIII), and superluminous supernovae (Leahy&Ouyed 

. .C^ , 2 AXPs, XTE J1810-197 and IE 1547.0-5408 which are ^^^^^ ^^^^^^^ ^J^^,^^ ^^^^ observed X- 

the only magnetars known to emit in the radio (Camilo et j j- j-- j ■ • j u j-- 

^ ray and radio properties, during quiescence and bursting 

al. 2006 . Both are demonstrably transient radio sources, phases 
having not been detected in previous surveys of adequate 
sensitivity. XTE J1810-197 is a transient AXfQ first de- 
tected when its X-ray flux increased ~ 100-fold compared 1.1. The X-ray emission 
to a quiescent level it maintained for at least 24 years 

(Ibrahim et al. 2004). Discovered with the Einstein X-ray In the pre-burst era, XTE J1810-197's ROSAT spectrum 

satellite in 1980, IE 1547.0-5408 was eventually identified showed a single blackbody (BB) with temperature T = 

as a magnetar candidate (Gelfand & Gaensler 2007) with 0.18 keV, an emitting area of - 520 km^((i/3.3 kpc)2, and 

spectral characteristics of an AXP. In this paper we look a luminosity of Lbb ~ 5.6 x 10^^ erg s~^{d/3.3 kpc)^ 

at these sources in the Quark- Nova context (hereafter QN; During its bursting phase, XTE J1810-197 showed a hot 

blackbody (T ^ 0.65 keV) with an exponential decay in X- 



* emaihouyedQphas.ucalgary.ca ray luminosity of ^ 280 days, as well as a warm blackbody 

^ In the sense that in quiescence their surface temperature (T ~ 0.3 keV) decaying at a rate of ~ 870 days (Gotthelf 
are as low as those of some ordinary young neutron stars & Halpern 2007). For the case of IE 1547.0 — 5408, after its 
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radio detection (Camilo et al. 2007a), an X-ray outburst 
was confirmed (Halpern et al. 2008) with a record high 
luminosity of ^ 1.7 x 10^^ (d/9.9 kpc)^ erg s^^ and with 
a total outburst energy of 10^^ erg < i?b < lO''^ erg. 

1.2. The radio emission 

For XTE J1810— 197, the radio emission began within 1 
yr of its only known X-ray outburst (Camilo et al. 2006 
and references therein). At its observed peak more than 
3 yr after the X-ray outburst, the radio flux density was 
more than 50 times the pre-burst upper limit. The X- 
ray flux has since returned to its quiescent level nearly 4 
yrs after the burst. IE 1547.0—5408, although not as well 
sampled as XTE J1810— 197, exhibits similar variations in 
flux density and was reported with a factor of 16 times the 
pre-burst upper limit (Camilo et al 2007a). 

Trends in radio emission between the 2 sources can be 
summarized as follows: 

— Both are very highly linearly polarized showing a flat 
spectrum over a wide range of frequencies. Their strik- 
ing spectra (i.e. spectral index > —0.5) clearly distin- 
guishable from ordinary radio pulsars (with a spectral 
index ~ -1.6; Camilo et al. 2007a&b). 

— At their peak, both magnetars are very lumi- 
nous in radio with luminosity at 1.4 GHz L1.4 > 
lOOmJy (d/kpc)^, which is larger than the L1.4 < 
2m Jy (d/kpc) of most any ordinary young pulsar 
(e.g. Camilo et al. 2006). For XTE J1810-197, its as- 
sumed isotropic radio luminosity up to 42 GHz is about 
2 X 10^0 erg s"! (Camilo et al. 2006). 

— Both have variable pulse profiles (exhibiting sudden 
changes in radio pulse shape) and radio flux densi- 
ties. The flux changes at all frequencies. At a given 
frequency there is no stable average pulse profile. 
Different pulse components change in relative intensity 
and new components sometimes appear. Sub-pulses 
with typical width approximately < 10 ms are ob- 
served (Camilo et al. 2007a&b). 

— For XTE J1810— 197, the torque was decreasing, at a 
time when the star was returning to quiescence years 
after the large outburst. As the torque decreased, so 
did the radio flux (Camilo et al. 2007b). 

— In IE 1547.0—5408, in contrast, the torque has been 
increasing, at a time when the X-ray flux has been 
gradually decreasing (Camilo et al. 2007a). 

In this paper we extend our existing quark star model 
for AXPs to account for the observed behavior of these two 
transients. This paper is structured as follows: Section 2 
gives the basic elements of the model. Section 3 describes 
the quiescent phase. Section 4, the bursting phase. Section 
5, the radio emission. Model predictions are highlighted in 
Section 6 before we conclude in Section 7. 



2. Basic components of the model 

The Quark- Nova is an explosive transition from a neutron 
star (NS) to quark star (QS) (Ouyed, Dey, & Dey 2002; 
Keranen&Ouyed 2003). The result is a partial ejection of 
the NS crust (Keranen, Ouyed, & Jaikumar 2005) that 
leads to two possible types of debris surrounding the com- 
pact remnant (i.e. the QS) depending on the QS's birth 
period. In OLNI, we showed that if the QS is born slowly 
rotating, then the debris formed from the QN ejecta will 
be in co-rotation, which we argue is responsible for SGRs. 
In OLNII, we showed that for QS born with millisecond 
periods, the debris evolves into a Keplerian ring with ap- 
plications to AXPs. Furthermore, RRATs are the result of 
late evolution of the Keplerian ring in our model (OLNIII). 
The interested reader is referred to these papers for more 
details. Below we give an overview of the salient features 
of the model in the case of a Keplerian ring before we 
apply our model to transient AXPs. 

2.1. Tlie quarl< star: magnetically aligned rotator 

The QN compact remnant is a quark star in the Color- 
Flavor-Locked (CFL) phase, which due to it's rigorously 
electric neutrality (Rajagopal & Wilczek 2001) possesses 
no crust. Owing to the superconductivity of the CFL state, 
the star's interior contains a lattice of vortices that confine 
the magnetic field (Ouyed et al. 2004). This interior con- 
figuration consequently forces the exterior field to align 
with the rotation axis (Ouyed et al. 2006; Niebergal et al. 
2006); this is illustrated in Figure [TJ 

2.2. The Keplerian ring 

Here we are concerned with compact remnants born with 
millisecond periods. As shown in OLNII, the quark star is 
surrounded by Keplerian debris we refer to as a ring. This 
Keplerian ring is described in detail in §2.1 and Appendix 
A in OLNII. Briefly, it is a high density (>~ 10^ g cm^^) 
ring rich in iron-group degenerate material (the ejected 
NS crust material) in Keplerian rotation around the QS. 
The ring is highly conducting, cool, and not threaded by 
the magnetic fleld. It is a (possibly amorphous) crystalline 
solid (like the outer crust of a neutron star). The ring ex- 
pands vertically and radially in time to a structure de- 
picted in Figure 2 with the ring inner radius at Ri^ ^15 
km, and an outer radius at i?out ~ 25 km (for these fidu- 
cial values the total area of the ring is Aj-i^g ^ 3000 km^ ) . 

The ring thickness in the z-direction can be shown to 
be i?ring ^ 2.68 kmp^^fg gi?ij^^ where /Oring,9 is the ring's 
density in units of 10^ gm cm~'^. In reality, on finer scales 
than depicted in Figure 2, the ring is subject to tidal frac- 
ture and is made up of many cylinders we refer to as 
"walls". The width of each cylinder is set by Keplerian 
shear resulting in meter size pieces. The mass and width 
of a wall are given in OLNII and are mw ^ lO~^"Af0 and 

3/2 

(5rw ~ 400 cm R-^ -^g, respectively. 
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Fig. 1. Illustrated here are the stages involved in the transition from a neutron star (NS) (inclined rotator) to a 
quark star (QS) (aligned rotator) in the Quark-Nova scenario. The collapse of the quark core induces an explosive 
NS-to-QS transition ejecting iron-rich degenerate crust material. The QS enters a superconductive phase confining the 
interior field to vortices, forcing the exterior field to align with the rotation axis (see Ouyed et al. 2004 for more details 
and Ouyed et al. 2006 for the related simulations). The iron-rich degenerate ejecta evolves into a Keplerian ring (at 
about 15-30 km from the star) surrounded, but not penetrated, by the dipole field (a co-rotating ejecta is also possible 
depending on the NS period; see OLNI). 



The magnetic field penetrates a conductor of thick- 
ness Sr on timescale tb ^ (Ana/c^) x (Sr)'^ where a is 
the conductivity (e.g. §4.1 in OLNII). The QS dipole field 
will penetrate the ring via its innermost wall. The wall is 
penetrated radially because Jr^ << ffiing- This occurs on 
timescales of a few hundred years (see eq.(17) in OLNII). 
As the QS magnetic field penetrates the innermost wall, 
magnetic torques (due to induced B^) slow down the wall, 
so it can accrete (see 2]). 



2.3. The ring atmosphere 

The ring's density decreases with height above the equa- 
torial plane. At any given temperature the ring's density 
below which the ring's matter becomes non-degenerate 
is found by equating the ring temperature to its Fermi 
temperature; this defines the ring's vertical atmosphere. 
The atmosphere is characterized by its base density, scale 
height in the z-direction (vertical to the orbital plane), 
scale height in the radial direction (i..e. on the outer ring 
edge, -Rout; see appendix B in OLNII), and thermal speed. 



P&tm 



460 gm cm-3 
20.4 cm 



(1) 



M3.3 



3.7 cm 



TkcvR 



^out,15 



M3.3 



9.4 X 10^ cm s~i 

^3.3 



where Tkov is the atmosphere temperature in keV, -R15 
is the radial position in units of 15 km, and /13.3 is the 
mean molecular weight in units of 3.3 which represents a 
partially- ionized iron- rich atmosphere (see §3.2 in OLNII). 
The ring atmosphere's scale height in the 2-direction is to 
be differentiated from that of the ring itself -ffring which 
is of the order of a few kilometers. 

The ring's atmosphere is penetrated by the magnetic 

field on timescales of ~ 12 days T^eV,o.i^n,i5kmf^llL,3.3- 
This can be derived from tb (as defined §2.2) using 
6r ~ -ffftm ^^'^ °^ l/cs where Cg oc \/ [ijT is the atmo- 
sphere's sound speed. The ring's temperature is in units 
of 0.1 keV representative of its equilibrium temperature 
during the quiescent phase; see eq(16) in OLNII. Since 
S^/Stt >> jOatm.^K, the atmosphere is forced to co-rotate 
with the field inhibiting accretion onto the star during 
the quiescent phase. On timescales of a hundred years, 
outbursts are triggered by wall accretion. These move de- 
generate ring material into the atmosphere (see §4 below) 
faster than the magnetic field lines penetrate vertically 
into the degenerate ring. The ring's penetration timescale 
is of the order of tens of millions of years because of the 
((5r)^ dependence, so that the ring is not penetrated from 
above. 

In our model, as we describe in more detail in §6.1 
in this paper, transient AXPs do not accrete during their 
quiescent phase while normal AXPs do accrete continu- 
ously from the ring's outer edge, where B^/Stt < Patm. 
For now we concern ourselves with transient AXPs whose 
quiescent phase is dominated by emission from vortex ex- 
pulsion as described next. 
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Fig. 2. The upper panel illustrates the ring structure during outburst. Fiducial values of the ring inner radius 
(15 km) and outer radius (25 km) are shown; the ring's vertical scale height is a few kilometers (see §2.2; the figure 
is not to scale). The ring is surrounded, but not penetrated, by the dipole field. The Bohm Front - "bump" made 
of non-degenerate Keplerian material - propagates outward at a speed Vbf heating up the system from Tq to Tb 
and inducing accretion onto the star. The quark star (an aligned rotator) is also shown with the accretion hot spot 
shrinking towards the pole at a speed ^hs a-s a response to the outward motion of the accreting Bohm front. The 
lower panel illustrates motion and fiow dynamics in the vicinity of the Bohm Front (BF). The solid line represents 
the transition to non-degenerate densities in the ring, pnd, while dsF is the radial width of the BF. As non-degenerate 
Keplerian material is expelled from the slow point, more Keplerian material is supplied from underneath (see Appendix 
[B|) . Accretion shuts off once the magnetic field penetrates the non-degenerate Keplerian layer re-enforcing co-rotation 
(behind the BF). 
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3. The quiescent phase in our model 

There are two critical radii in our model during the quies- 
cent phase, the inner radius i?in, and the outer ring radius 
i?out- In most cases these radii will be expressed in units 
of 15 km or 25 km thus assigned a subscript 15 or 25. The 
other 2 parameters related to the geometry of the ring 
are the ring's solid angle divided by 47r at i?in and i?out 

namely, /in,n = -ffringan/-Rin and /out.O — ^^ring,out/-Rout- 

General relativistic (GR) effects are included in the factors 
/gr,qs = \/l — Rsch./Rqs, Jgrm = ^ RscU./ Rin, 

/cR.out = \/l — ^Sch./^out while /cR.ring — 0.5(/GR,in + 

/cR.out) with i?sch. being the star's Schwarszhild radius. 
Unless otherwise specified, quantities such as luminosity, 
temperature, and area are local values. Values at infinity 
are obtained by using the relevant GR factors. Finally, the 
ring area which includes top, bottom and inner surfaces 

is Aout = 27r(i?2^t - RfJ + 47ri?ini?ring,in = 2^^Rl^^Yout 



lout = ^-^[.Rlut 

where y:„t = 1 



ing,in — 

{Rin/ Rout)^ + 2(i?in/ Rout)^ fin,n- 



3.1. The 2 blackbodies in quiescence 

As discussed in OLNI and OLNII, during the quiescent 
phase we have 2 blackbodies, one from the emission due 
to magnetic reconnection following vortex expulsion (the 
emission occurs just outside the star's surface). The re- 
sulting luminosity is (see §5 in OLNI) 



ivortox 2 X 10^* erg s ^x.o.i-P-n 



with a corresponding temperature. 



0.2 keV 77x(o. 



P 



1/2 
-11 



-^QS,10 



(2) 



(3) 



where the period derivative is in units of 10~^^ s s~^, the 
star's radius in units of 10 km, and rjx is the efficiency pa- 
rameter inherent in the conversion from magnetic energy 
to radiation in units of 0.1. 

The second BB results from reprocessing by the ring 
of the first BB's X-ray emission, /qr Qg/in.oivortcx = 
/gr ring^ringC''^*- The ring-atmosphcrc system's temper- 
ature during quiescence is then, since /gr,qs — /or, ring, 



To 



pl/2 



p3/8 ^1/4 
^out,25^out 



(4) 



where the subscript "q" stands for quiescent in contrast to 
the values during the bursting phase denoted by subscript 
"b" . These two blackbodies are generic emission compo- 
nents to reprocess the emission from the central object, ir- 
respective of its ultimate origin. As such, they should also 
be expected in normal AXPs/SGRs, where the continu- 
ous emission is dominated by the constant accretion from 
the ring edge rather than by vortex annihilation. Since 
P oc in our model (Niebergal et al. 2006), both BB 

temperatures during the quiescent phase evolve in time as 



The resulting temperature ratio in our model, using a 
fiducial value l^ut ~ 0.64, is 



1.7- 



3/8 
out, 25 



1/2 
QS.IO 



(5) 



Interestingly a similar correlation between the hot 
and cool BB for AXPs and SGRs has been ob- 
served (Nakagawa et al. 2009; see their figure 4), with 
Tbb.h/Tbb.c ~ 2.5. This we suggest as evidence for re- 
processing in the system. The list of objects studied by 
Nagakawa et al. (2009) does not include transient AXPs 
in quiescence. These seem to be modeled by a single ther- 
mal component (e.g. Gotthelf & Halpern 2005; Perna & 
Gotthelf 2008; Bernardini ct al. 2009), although because 
of the low fiuxes, presence of the cooler component cannot 
be ruled out. 



4. Bursting phase in our model 

4.1. Consequences of wall penetration and accretion 

The ring remains quiescent until magnetic penetration of 
the inner edge of the ring (wall) . The magnetic field pene- 
tration and subsequent accretion events are very sporadic 
- they last for hours and occur roughly once every hundred 
years (see section 4 in OLNII) . The consequences of wall 
accretion are as follows: 

— Change in ring's mean molecular weight: 

The energy released by the wall accretion is sufficient 
to dissociate a significant mass of iron nuclei in the ring 
into light nuclei (nuclei with Z ~ 13). The state of the 
ring depends on the Coulomb parameter which gives 

1 /3 

a sohdification temperature of Tg ~ 9.5pg keV (see 
§3.2 in Ouyed & Leahy 2009). The ring mean density 
(eqn A.8 in OLNII ) is - lO^t^/^ ~ 3 x 10^ g cm'^, 
with age t in yrs. This gives Tg c± 1.4 keV. However 
during the wall accretion event, the ring temperature 
rises to ^ 4-5 keV (see eqn C.2 in OLNII). This com- 
pletely melts the ring, which allows light elements to 
rise buoyantly. Then the ring re-solidifies on a few hour 
timescale (eqn 22 of OLNII). This creates a two layer 
system as depicted in Figure 2. This process reduces 
the molecular weight of the atmosphere from /iq 3.3 
to /j,b ~ 2.1-2.5 (we adopt an average /ib ^ 2.3). 
As estimated in §4.4 in OLNII, for a typical wall 
mass Tow ^ lO~^°Af0, the number of dissociations 
following irradiation from wall accretion is 7Vi3 
10'"'Cw,o.ooi?7o.iWw,-io where Cw,o.ooi is the dissocia- 
tion efficiency in units of 0.001, and ryo.i is the wall 
accretion efficiency in units of 0.1. The depth of the 



Z ^ 13 layer can be estimated to be H13 ^ 10~'^iJring; 
that is, of the order of a few meters. Subsequent deple- 
tion of the Z ~ 13 nuclei by accretion leads to a return 
to an iron-rich atmosphere (i.e. fj,q 3.3) as shown in 
OLNII (see also Appendix IB. 21 in this paper) . 
The Bohm front (BF): 
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Fig. 3. Upper panels: Model luminosity (i/lO"^^), temperature and area (curves) from one side of the ring as 
observed at infinity compared to observations (dots) of the warm component of XTE J1810— 197. The dotted line is 
the contribution inward of the BF, the dashed is from outward of the BF, and the solid is the combined contribution 
(see text); Lower panels: Model luminosity (L/IO^*), temperature and area (curves) from one pole as observed at 
infinity compared to observations (dots) of the hot component of XTE J1810— 197. In the right-most panel the dotted 
line traces the shrinkage of the HS as the BF moves outward, the dashed line shows the late contribution from accretion 
induced by heating of the back side of the ring, and the solid is the combined contribution. 



The sudden reheating of the inner ring region triggers 
heat propagation outward along the ring's atmosphere. 
Heat propagation normal to the magnetic field can oc- 
cur by either classical or Bohm diffusion. 
The classical diffusion coefficient is D± = 
2nr]j_kBTc^/B'^ where ri± = 1.44 x 10-8ZT-^/'^ln{A) 
is the transverse Spitzer resistivity with ln(A) the 
Coulomb logarithm (e.g. eq. (5.71) in Chen 1984 for 
S. I. units). The Bohm diffusion coefficient is given 
by ckBT/{16eB) (see eq(5.111) in Chen 1984 for S. I. 
units). The ratio of classical to Bohm coefficients for 
B = 10^3 G, T = 0.2 keV and density of 10 g cm^^ 
is of the order of 10^® implying that Bohm diffusion 
dominates over classical diffusion for magnetic field 
strength and temperatures involved here. Thus the 
heat diffuses outward according to Bohm diffusion, 
introducing a critical radius in our model (during the 
bursting phase), Rbf = Rin + Ar, where Ar is the 



distance that the Bohm Front (BF) has travelled from 
Rin (see Appendix A for details). 

— Atmosphere feeding: 

The boundary between the non-degenerate atmo- 
sphere and degenerate ring material, at the BF, moves 
downwards into higher density layers since the BF 
heats up the interface to higher temperatures (re- 
call that Patm ~ 2^keV' eq.(IT|)). Thus newly non- 
degenerate Keplerian ring material is unveiled between 
the CO- rotating atmosphere and the interface (Fig. 2). 
We thus have two mechanisms that feed the atmo- 
sphere with new non-degenerate material. These two 
contributions are seen when estimating the surface 
density of the atmosphere S Patm^^atm oc T^^^/ n; 
the lower /i from the buoyancy and the higher T from 
the heat front will increase the mass up by a factor of 
- 10. 

— Atmosphere ejection and accretion: 
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Tlic huge shear between the co-rotating, magnetized 
atmosphere and the underlying non-degenerate keple- 
rian atmosphere leads to Kelvin-Helmholtz instability 
that helps load the keplerian material onto the mag- 
netic field lines attached to the heat fronl[l. With sim- 
ple angular momentum arguments one can show that 
any keplerian particle attached to a co-rotating field 
line will slide along the field line moving radially out- 
ward. The magnetic field acts as ramp for the plasma 
particles to get rid of their excess angular momentum. 
This is the slingshot effect in magneto-hydrodynamic 
(MHD) jets (discussed in Appendix [B|) . In the quies- 
cent state there is a huge shear between the corotating 
field lines (containing the non-degenerate atmosphere) 
and the Keplerian degenerate ring material. However 
since the field lines do not penetrate the degenerate 
ring no MHD ejection (and thus no accretion) is feasi- 
ble. 

In summary, the K-H instability would load mass onto 
the field line attached to the Bohm front and get 
ejected centrifugally along the field lines as explained 
in Appendix IB] (see also below). The wind is then chan- 
neled onto the star with an accretion rate, mace, given 
by eauation (jB.2p creating a hopt spot (HS) on the sur- 
face of the star. 



4.2. The 2 blackbodies during burst 

Only magnetic field lines that are in the path of the heat 
front get loaded and accrete onto the star. The main con- 
sequence, as we show in details here, is that the HS on the 
star moves closer to the pole thus decreasing in area. 

The temperature of the inner ring is obtained by equat- 
ing heating from BF induced accretion onto the quark 
star with blackbody cooling by the inner ring surface 

(i? < Rbf)] /GR,Qs/BF,f2iacc. = /cR.ring^BFO-T* where 

-^^acc. — V^accC^ is the accretion luminosity, r] the ac- 
cretion efficiency, /bf,o = -faring, bf/-Rbf, and the area 
created by the Bohm front is Abf — 27r(_R|p — + 
47ri?inFringan = 27ri?|pYBF- Here Ybf = 1 - (i?i„/i?BF)^ + 
2(i?in/^BF)^/in.si- The resulting temperature is, 



Tb ~ 0.21 keV T^o.iC 



,sp,0.01 



RbF.15 



3/2 Yf 



BF 



(6) 



Here, Csp,o.oi in units of 0.01, is a factor related to mass- 
loading at the slow MHD point as determined in Appendix 
B. The approximation above is valid until late times when 
accretion ceases, so that Tb is no longer given by the 
accretion-cooling balance but rather settles to the qui- 
escent temperature Tq. The outer ring's temperature is 
given by the quiescent temperature, eq.(|4|), until the heat 
front reaches it. 



^ During quiescence the shear is between the solid degenerate 
Keplerian ring and the co-rotating atmosphere thus not prone 
to the instability. In contrast, during burst the shear is between 
one fluid (the new Keplerian atmosphere) and another fluid 
(the pre-existing co-rotating atmosphere). 



The accretion rate is then obtained by combining 
ea. (IB.2p with equation above to get 



1.6 X 10^5 g s-i %'.iCp,o.oi 



R 



11/2 
BFaS 



/K(rb) 



,,6y3 
P ^BF 



(7) 



The function /K(2b) is described in Appendix [Bl and acts 
as a "valve" that shuts-off accretion after the heating front 
has passed. There is a fixed amount of Keplerian material 
at any radius that can be fed into an MHD wind and 
thus accreted onto the star. This is because of re-enforced 
co-rotation once the magnetic field penetrates the non- 
degenerate Keplerian layer effectively shutting off accre- 
tion (see Figure 2). 

The natural connection between the BF outward ex- 
pansion and the resulting HS on the star is illustrated in 
Figure 2. The corresponding HS luminosity is 

Lace. - 1.4 X 10- erg s^ ^iC^o.oi^^^^^^ , (8) 



M ^BF 



with a HS temperature 



Ths ^ 0.58 keV ?7o.iCsp,o.oi 



,,3/2^3/4^1/4 
A* -^BF ^HS,120 



(9) 



where the HS area, Ahs is in units of 120 km^, is derived 
in the Appendix (see eq. IC.ll) . 

In summary, during the burst phase, the ring- 
atmosphere system consists of 2 BBs; the inner warmer 
part increasing in area at the Bohm diffusion rate and the 
outer cooler part decreasing in area. The hot BB in our 
model is provided by the HS which decreases in area as 
the BF propagates outwards. 



4.3. Temperature ratio during burst 

The temperature ratio in the bursting phase case is then: 



A 



BF 



a/4 



' /bf.^^hs ' 



2.8 



out, 25 



Qsao 



1/2 



(10) 



^ring/2, /bF,0 ^ /out, SI, 



where we approximated Abf 
and ^HS O.IAqs to get the numerical value. Our model 
predicts the temperature ratio during burst for a given 
object to be higher by a factor of ^ 1.5 than during quies- 
cence. We expect a scatter in comparing different objects 
mainly caused by variations in i?out- This is consistent 
with Figure 4 of Nakagawa et al. (2009). 

4.4. Application to XTE J1810-197 

The upper panels in Figure 3 compares the time evo- 
lution of the ring's luminosity, temperature and area in 
our model to the observed warm BB in XTE J1810-197. 
In our model the luminosity from the inner part of 
the ring is Lbf = AbfctT^ while the outer part gives 
Lout = AoutaT^ with Aout = 27r(i?2^t - Rip). The to- 
tal luminosity from the ring during burst is Lring,b = 



8 
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Lbf + ^out while the corresponding effective area and 
temperature were derived by weighting over luminosity, 

^ring.cff. = (^BF^BF + ^out-^out)/iring,b and Tringxff. = 
ring,b- 

The lower panels in Figure 3 compare the time evo- 
lution of the HS's luminosity, temperature and area in 
our model to the observed hot BB in XTE J1810--197. 
In this case, only one component comes into play, the ac- 
cretion luminosity Lace- However, as the BF gets closer 
to i?out the back side of the ring is heated resulting in 
additional accretion by the same mechanism (i.e. eating 
into Keplerian material in the radial direction this time) 
and additional area on the polar cap defined by the field 
lines connecting the back side of the ring to the star (see 
Figure 2). The resulting area is ^hs + ^cdgc and temper- 
ature (iacc./(CT(AHS +^cdgc)))^/''. 

The fits to XTE Jf 810-197 data were obtained for the 
following set of parameters: 

i?i„ = 23.5 km (11) 
i?out = 26.5 km , 

and by slightly adjusting the mass-load at the slow MHD 
point so that Cs-p — 0.02 (see discussion following eq. jlB.ip ). 
The other parameters were kept to their fiducial values 
including the star's parameters i?Qs = 10 km, Afqs — 
1.4Mq. The star's magnetic field is given by i?Qs = 

VSkPP ~ 3.3 X 10^^ G (with K = 8.8 X 10^^ s^^ 
as given in eq.(lO) in OLNIII). 

Assuming that XTE J1810— 197 has experienced a few 
bursting events, using equation (A. 7) in OLNII, the ring 
would have spread to no more than (Ar)t ~ 8 km if the 
system's temperature during quiescence remained on aver- 
age ~ 0.1 keV. This is consistent with (i?out — ^in) found 
from fits to the XTE J1810-197 data thus providing a 
self-consistency check on our model. This also confirms 
our overall findings in previous work (OLNI, OLNII and 
OLNIII) that the ring should be a few kilometers in width 
after a few hundred years. 

4.5. Application to IE 1547.0-5408 

In its bursting phase, this source was fitted with a T ^ 0.5 
keV hot BB with Lbb 1-3 x 10^^{d/9 kpc)^) erg s^^ and 
a corresponding area decreasing from 180 km^ in June 
2007 to 96 km2 in August 2007 (see Halpern et al. 2008). 
During the quiescent phase, Lbb 10'^'*(c?/9 kpc)^) erg 
s"\ Tbb ^ 0.4 keV and Abb ^ 36 km^ (see Table 1 in 
Halpern et al. 2008). 

This source has not been as well sampled as was XTE 
J1810— 197. Nevertheless, in our model, high BB tempera- 
tures during the quiescent phase are suggestive of a more 
compact ring which is closer to the star (see eq.(|4])). In 
our model, such a small and compact ring could have 
been a consequence of a smaller amount of crust material 
ejected during the QN (see eq.(2) in OLNII). For exam- 
ple, i?out = 15 km and i?in = 13 km, which implies l^ut 
Ybf ^ 0.117, inserted in eq.gl) gives Tq ^ 0.35 keV; we 



take P-u ^ 2s for this source. The ring's burst epoch 
temperature is then ~ 0.48 keV with a corresponding peak 
luminosity from eq.([H]) of Lace. ^ 4.8 x 10^^ erg s^^ for 
/ib = 2.4. This is close to the - 1.3 x 10^^{d/9 kpcf) 
erg s"""^ measured in June- July 2007 (the peak of the out- 
burst was not observed and could have been higher than 
this). The initial area of the HS is given by equation (IC.ip 
and is estimated to be ^ 80 km'^ using Rm = 13 km and 
^ 15 km. Finally, in the case of IE 1547.0-5408 there 
seems to be hints of an aligned rotator from its small X-ray 
pulsed fraction and its relatively broad radio pulse. In our 
model, the QN compact remnant (the QS) is born as an 
aligned rotator due to the vortex confinement of magnetic 
field (see Figure [T|). 

5. The radio emission in our model 

In this section we develop a scenario for radio emission. As 
discussed in §1.2, XTE J1810-197 and IE 1547.0-5408 
show unique characteristics in radio (including delay after 
X-ray outburst, a flat spectrum, and unusual spin-down 
behavior). In chronological order the sequence of events 
that leads to radio emission in our model is as follows: 

— Magnetic bubbles are generated at the ring during the 
X-ray outburst (see §5.1). 

— The propagation time of these bubbles from the site 
of production (the ring) to the I.e. defines the delay 
between radio and X-ray outbursts (see §5.2). 

— The unusual spin-down behavior is induced by torques 
on the magnetosphere from the co-rotating bubbles, 
during their outward propagation (see §5.5). 

— Relativistic collisionless reconnection at the I.e. de- 
stroys the bubbles. 

— Magnetic energy released by the bubbles is eventually 
radiated by particles accelerated by reconnection at 
the I.e., mostly at radio frequencies (§5.3). We suggest 
that the flat radio spectrum is naturally associated to 
this mechanism (see §5.4 below). 

5.1. Magnetic reconnection and bubble generation 

After penetration, the poloidal magnetic field (Bp) lines 
inside the wall (inner ring) are dragged by the Keplerian 
shear generating a toroidal magnetic field, B^. Continuous 
reconnection events during the winding of the field lines 
lead to X-point generation and emergence of closed mag- 
netic bubbles; more concisely these magnetically confined 
plasma bubbles (plasmoids) result from the Keplerian 
shear. 

Generation of bubbles is a common feature of threaded 
disks as seen in many simulations (e.g. Romanova et al. 
1998; Yelenina et al. 2006). These simulations show the 
magnetic bubbles to acquire enough speed to escape grav- 
ity and expand freely outward. Magnetic loops and bub- 
bles are also common in the sun (referred to as plasmoids) 
and are also found to gain enough energy from the recon- 
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nection events to escape the system and expand freely (e.g. 
Wagner 1984; Tamano 1991). 

One key difference between these cases and our model 
is the fact that in our model, during quiescence, the de- 
generate ring is not threaded by the magnetic field. The 
bubbles form only during the bursting phase once the mag- 
netic field has penetrated the inner ring and is sheared. 
The B field attached to the broken inner pieces of the wall 
is wound up by the keplerian angular velocity {'^ 6000 
rad/s) of the piece vs. the corotation angular velocity (~ 1 
rad/s) of the footpoint of the fieldline on the quark star. 
This rapid field line winding results in reconnection and 
magnetic loop (bubble) formation. The B field penetra- 
tion (and subsequent wall accretion with concurrent bub- 
ble generation) events are very sporadic and unique to our 
model - they last for about an hour and occur once every 
century (see §4 in OLNII). 

A rough estimate of the number of bubbles that can 
be generated is 

iVBub. ^ ~ 8.3 X lO*'^^,^ , (12) 

I'Bub. 

where the total reconnection volume is Viv ~ 
47ri?ini?ring and Veub. ^ (47''/3)5r^. The mass 

of a given bubble is thus rriBub. ^ (47r/3)(5r^Patm.,w 
where the plasma confined by the bubbles has a den- 
sity at birth given by the atmosphere density Patm.,w 
The magnetic energy stored in each bubble is then 
EBuh. ~ (47r/3)(5r3 (52/87r) with ~ Bp,i„ with Bp,i„ = 

^Qs(^Qs/^in)^- 

5.2. Delay between the X-ray and radio emission 

Simulations of bubbles in threaded disks show that mag- 
netic loops and bubbles propagate outward because of 
the energy they gain from reconnection and because of 
the dipolar B-field gradient from the central star, simi- 
lar to the way solar field disturbances propagate outward 
above the solar photosphere. These bubbles propagate at 
the sound speed of non-dissipative magnetic field distur- 
bances, ie. at the Alfven speed. Illustrated in Figure 4 
in this paper, is a rudimentary representation of what is 
seen in these simulations; an outward migration process 
through the dipole field involving reconnection. This is a 
simple and only qualitative model for bubble propagation, 
which is a propagating geometric disturbance in the mag- 
netic field/plasma. Bubble propagation results in different 
field lines and different plasma making up the bubble as 
it moves. As seen from figure 4, a magnetic field line in 
front of the bubble is only temporarily part of the bub- 
ble, and leaves the bubble after the bubble has propagated 
past it. The reconnection on the front side of the bubble 
results in surface pressure similar to the static B'^ /8tt, so 
the dynamics of reconnection during bubble migration are 
important to how the bubble evolves (both its speed and 
change in volume) as it migrates. 

The Alfven velocity depends on the magnetic field ge- 
ometry inside the light cylinder (which we take to be 



nearly dipolar) and on the ambient density which we take 
to scale as cx (iiin/r)". If there is a steady wind inside 
the magnetosphere, then if the wind is spherically sym- 
metric the density we expect would decline as or if it 
is equatorially confined the density would decline as r~^; 
that is 1 < a < 2. Although the exact distribution is 
not at all known (e.g. Michel 1969; Spitkovsky 2008; see 
more discussion on this aspect of our model in l6.5|) we will 
adopt a = 1 as we expect the bubbles to remain along the 
equator since they are Keplerian at birth. 

The bubbles move outwards at a rate given by 
dr/dt = VA where va = c{Ro/r)^~°'^^ so that Rq ^ 

/ 5/4 Nl/(3-a/2) 

45 km ( BQS44/Tati„^kcV ) the radms beyond 

which the Alfven speed becomes sub-relativistic; the star's 
magnetic field i?Qs,i4 is in units of 10^^ G. 

The time it takes the magnetic bubbles to reach the 
light cylinder (Ic) at R\c = c/f2 (at which point they start 
dissipating) is found from Jjj^'" dr/vA, 

where we made use of the definition of i?o and Bqs oc 
\/ Pp in our model. For example for a = 1, we get 

p32n5/4 

V-op- ^ 4.1 yrs r*"^''"'^ , (14) 

-f-ii 

where the period is given in units of 5 seconds; recall that 

Sqg = 3.6 X 10^^ G \J P5P-11. As they slow down while 
flowing towards the Ic the bubbles will pile-up. This im- 
plies a delay between the X-ray outburst following wall 
accretion and the peak of radio emission in our model. 
Thus, following the peak of the X-ray outburst for XTE 
J1810-197 (with P ~ 5.54 s) and IE 1547.0-5408 (with 
P ~ 2 s), we expect a corresponding radio peak emis- 
sion to occur with a delay of ^ 320 days and ^15 days 
respectively (assuming Tatm. ^ 0.3 keV and P-n ~ 2). 

5.3. Duration and luminosity 

The radio duration has two components in our model: 
(i) The radiative lifetime of the electrons, ^cooi; (n) The 
time difference between when the first bubble arrives and 
the last bubble arrives at the I.e., ^radio = Aipj-op.. As 
we show at the end of this section, propagation delays 
between bubbles is the dominant component. Below we 
focus on case (ii). 

Since all bubbles are produced within a few hours 
(during the X-ray burst) and at the same location, the 
radio duration is determined by variations in propaga- 
tions times of the bubbles to the I.e.. The magnetic re- 
connection events lead to variations in B which translate 
to fluctuations in va and thus propagation time (con- 
sequently arrival time at the I.e.) resulting in iradio = 
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Fig. 4. Keplerian shear during wall accretion leads to X-point reconnection which in turn causes the generation of 
closed magnetic loops, or "bubbles" (see §5.1 in text). In our model, these magnetic bubbles are responsible for the 
radio emission seen in some AXPs. The top panel illustrates the outward migration of a bubble through the surrounding 
magnetic field. The bubble's magnetic polarity in front as it moves outward is opposite to that of the surrounding 
dipole. This leads to reconnection and the motion of the bubble towards the light cylinder. The lower diagram shows 
the overall pattern of bubble migration (i.e. magnetic buoyancy) and expansion that leads to pile up near the light 
cylinder. The bubbles dissipate as they cross the light cylinder leading to the radio emission (see text for details). 



At 



prop. 



it 



prop., max ^prop. ,min 



^prop 



(AB/B) - (B,„ax - S,ni„)/S. Or, 

AB 

^radit 



i) = iprop X i^B/B) where 



iprop. X ^ 



(15) 



For a homologous expansion during the outward prop- 
agation, a given bubble increases in size as SrBuh. ~ 
^^w(^/^in)- Combined with flux conservation, B^nh. oc 



5r 



Bub. 



and Veub. oc Sr%^^^, it yields E'sub. 



VBuh.B'^uh /(8^) oc r ^. The radio luminosity is then 
given as, for a ~ 1, 



L 



radio 



TT^Bub. 
mag. , I.e. 

^ radio 



(16) 



B 



4.4 X 10^® erg s x 7?r x - — x 

p3/2 n6 „l/6 
^-ll^QS,10^'ring,9 

-^5 -f^m,15-tatm.,kcV 

where E^H,^ = E^^.^^ x (R^jRi.,.) ~ {%^V^) x 
(i?in/i?i.c.) is the total magnetic energy stored in the bub- 
bles by the time they reach the light cylinder. The ratio 



(i?in/i?i.c.) is the dilution factor induced by the homolo- 
gous expansion of the bubbles as they propagate from Rm 
to i?i.c.; recall that Bin — Bqs{Rqs/ Rin)^ ■ In the equa- 
tion above, 77R is the efficiency of conversion of magnetic 
energy to radio emission via reconnection (e.g. Ouyed et 
al. 2006) at the I.e.. 

For XTE J1810-197, Camilo et al. (2007c) report radio 
emission lasting at least 272 days with initial fading from 
5 mjy to 2mJy in the first few weeks (see their figure 
2). In our model, tj-adio ^ 272 days gives AB/B ~ 0.85. 
Then using i?in — 23.5 km (from the X-ray fits in §4.4), 
iradio — '7R X 3 X 10^° erg s~^. This matches the observed 
radio luminosity of ~ 2 x 10^° erg s~^ for rj^ ~ 0.7 which 
is suggestive of a high dissipation/reconnection efficiency 
of the bubbles at the I.e.. 



For IE 1547-5408, radio emission was observed 
(Camilo et al. 2008) from June through August 2007 (~ 90 
days). It was not detected Jan 22, 2009 within 18 hours of 
the first reports of renewed X-ray (SCR) bursting activity 
from it (Camilo et al. 2009). In our model, this is sugges- 
tive of AB/B - 6.0. Then using i?i„ ~ 13 km (from the 
X-ray fits in §4.5), iradio — x 2 x lO'^^ erg s~^. This 
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matches the observed radio luminosity of ^ 2 x 10'^'^ erg 
s-i for 77R ~ 1%. 

The fits above indicate an order of magnitude differ- 
ence in efficiency, 77R, between the two sources. However, 
a smaller a (i.e. a slower density decrease inside the I.e.) 
for IE 1547-5408 would lead to a longer propagation time 
and lower radio luminosity, implying a higher reconnec- 
tion efficiency. Alternatively, a more compact quark star 
for IE 1547-5408 would also imply a higher reconnection 
efficiency. 

As a given bubble reaches the I.e. it is sheared and 
spread-out around the equator of the I.e.. The result- 
ing bubble's magnetic field is given by BB^b. _ 



ryBub. 



,Biib. 
I.e. 



60 G 



X 27ri?i.c.. This yields 



pl/2 n.3 



QS,10 



p3/2pl/2 
^5 -^in,15 



(17) 



After shearing, reconnection occurs accelerating electrons 
(see tj5.4l) . The electrons emit synchrotron radiation with 
a critical frequency of Vc = 0.42GHz7j'Qi?|^" q where 710 is 
the electron Lorentz factor in units of 10 following recon- 
nection and the sheared bubble magnetic field is in units 
of Gauss. We get 



25 GHz 



pl/2 n3 
^-ll^C 



QS.IO 



p3/2nl/2 
-^5 ^in,15 



(18) 



For XTE J1810-159 the above gives ^ 24 GHz while it 
is Vc 150 GHz for IE 1547-5408. Observations at higher 
frequencies than currently available could constrain our 
model. 

The above calculations neglect the radiative lifetime 
of the electrons. However in cases where the bubbles 
reach the I.e. simultaneously (i.e. with negligibly short 
time delays), the radio duration will be dominated by 
the radiative cooling lifetime of the electrons, tcooi ^ 
0.25 day/(7ioi?sij gg) (Lang 1999). For our estimated val- 
ues of 7 and Bsh_, , the propagation delay is much longer 
than the electron radiative lifetime. E.g. for IE 1547- 
5408, the radiative lifetime only becomes important, with 
7 — 10, for Bs\i. < 3 G, which would also result in the 
critical synchrotron below the observed radio frequen- 
cies. 

5.4. The flat spectrum 

As they cross the Ic, the bubbles will dissipate by brak- 
ing up into smaller bubbles while driving complex, non- 
linear reconnection events. Particle in Cell simulations of 
X-point reconnection events in a pair plasma have shown 
that highly variable radio emission, with an extremely fiat 
spectrum (s —0.17) is a natural outcome with pairs ac- 
celerated to Lorentz factors of up to ^ 70 (e.g. Figure 2 
in Jaroschek et al. 2004). The flat spectrum is observed 
up to the cutoff frequency at ^ 100 GHz (e.g. Figure 4 
in Jaroschek et al. 2004). These simulations show that 
the highly dynamic non-linear evolution of thin current 



sheets serve as the fundamental plasma scenario to ob- 
tain a flat synchrotron power spectra in pair-dominated 
environments. The emitted synchrotron emission is highly 
linearly polarized, with spectral polarization varying from 
> 50% around 1 GHz and up to 80% or higher at higher 
frequencies. 

If pair generation regions exist in the vicinity of the Ic 
then, when combined with magnetic reconnection events 
from bubble dissipation, it would offer a natural explana- 
tion for the observed flat spectrum. This could also help 
account for the prevalence of emission at a particular rota- 
tion phase (see fig. 2 in Camilo et al. 2006) which implies 
that the co-rotating bubbles must preferentially dissipate 
at a particular co-rotating longitude. It might also be the 
case that there exist a small non-uniformity at the light 
cylinder, caused by feedback from pairs produced by bub- 
ble dissipation mechanism itself. 

5.5. Torques during radio emission 

As each bubble expands toward the light cylinder it ex- 
erts a torque Tsub. = x diBuh./dt = fl x {2mBuh.TVA + 
"rriBub.r^) where /Bub. = mBuh.r^ is the moment of inertia 
of a bubble at a radius r from the star; here rheub. is the 
rate of change of the bubble mass. For riiBub. = 0, the 
total torque exerted, t — A^Bub.TBub. , is then 



r ^ 8.8 X 10^2 dyn cm^-^^tf^' 



''days 



(19) 



The corresponding frequency derivative, v = —t/{2t:Iqs) 
^QS 



with /qs ~ lO''^ g cm^ for the star, is 



-1.4 X 10"^^ s"2 



-"•in, 15 ,-3/7 
P5IQSA5 "^^^ 



(20) 



This gives ~ —1.9 x 10 



and 



-10 



-13 „-2 



for 



XTE J1810-197 and IE 1547.0-5408, respectively, which 
is of the right order of magnitude compared to what is 
observed. 

In our model, the torque decreases in time (as the X- 
ray decreases) in accordance with what has been observed 
in the case of XTE J1810-197 (Camilo et al. 2007c). 
However, for IE 1547.0—5408, the torque increased in time 
as the X-ray flux decreased. To explain the case of IE 
1547.0—5408 we recall that the first term {fl x 2mBuh.'i'VA) 
in the torque equation gives a decreasing torque with ra- 
dius while the second one {fl x rh-Bub.''^) gives an increas- 
ing torque with radius. Since the radius increases in time 
this is equivalent to the torque changing accordingly in 
time. In order for the torque to increase with time then 
rnBuh.T^ > 2mBub.™A which at a distance half way to the 
I.e. implies mBub./"^Bub. > 10^^/P''/^. Thus a slight in- 
crease in the bubble's mass during propagation can give 
an increasing torque with time. 
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6. Discussion 

6.1. Transient versus normal AXPs during quiescence 
in our model 

The two sources studied here are termed transients in the 
sense that in quiescence their measured temperatures are 
as low as those of some ordinary young neutron stars. 
Only during their bursting phase that their bolometric lu- 
minosity becomes comparable to that of a typical AXP in 
quiescence (i.e. a bolometric luminosity of ~ 10"^^ erg s^^; 
see Figure 2 in OLNIII). From approximate outburst time, 
the estimated initial bolometric luminosities for the two 
transients are ~ 2 x 10'^^ (i| 3 erg s~^. 

In our model, to understand the difference between 
transient and normal AXPs during quiescence, one should 
note that the magnetic field would impose co-rotation of 
the ring's atmosphere (thus no accretion onto the star) as 
long as the magnetic energy density anywhere along the 
ring's atmosphere exceeds the Keplerian energy density. 
The condition is 

Thus during the evolution of the source as long as the 
condition above is satisfied there will be no accretion and 
the object's emission is dominated by vortex expulsion. 

The star's magnetic field decays at a rate (see §3 in 
Niebergal, Ouyed, & Leahy 2006) 

Bs = Bo[l + -j (22) 

where T = 840 s is the characteristic age due to vor- 
tex expulsion using our fiducial value of Bqs,q = 10^^ G 
and Pqs.o = 1 nis for the star's surface magnetic field 
and period at birth, respectively. Simultaneously, the ring 
spreads viscously outward (i.e. radially) at a rate given by 
eq.(A.7) in OLNII 

(Ar)ri„g ^ 7.8 kmr^'/v^yr' ' (23) 

where time is in years. The ring's radius at birth is given 
by equation (2) in OLNII and is i?out,o — 15 km/m^y 
where m_7 is the ring's mass is in units of IO^^Mq so that 
it spreads from i?out,o to i?out = -Rout,o + (Ar)i.ing; clearly 
for sources thousands of years old (as is the case here) the 
ring's outer radius is simply given by i?out ^ (Ar)i-ing. The 
condition for no accretion, C < 1, is then equivalent to 

t < 1300 yrs x Tjf/'' x S^^f . (24) 

To get equation (24) we start with condition (21) and 
i?out = Arring (i.e. eq. 23) and replaced Bg from eq. (22) 
with t » T] The QS mass and radius are kept to their 
fiducial values. This implies that for Bq = 10^^ G and 
0.2 < T(keV) < 0.4, the transient age is 3 x lO'' < t(yrs) < 
4 X 10'^. Transient AXPs, we speculate, should be younger 
than accretion-dominated sources (such as 1E2259-I-586 
and 4U0142-H61) in our model. 



The magnetic field decay from vortex expulsion im- 
plies that magnetic field lines defiate radially inwardfH, 
so we expect that violation of co-rotation to first occur 
at the ring's outer edge. Also, at the outer edge of the 
disk, there is magnetic field of opposite B^ being pro- 
duced when matter drags field lines. This is an unstable 
situation since reconnection can occur between adjacent 
field lines of opposite polarity leading to accretion along 
the reconnecting field lines. Anywhere else on the surface 
of the disk, whenever C > Ij the atmosphere winds up the 
field lines without reconnecting them until the magnitude 
of the magnetic field is strong enough to re-inforce coro- 
tation, ie. the magnetic field stiffens, leading to a stable 
balance. 

In summary, a transient AXP quiescent phase is dom- 
inated by X-ray emission from vortex expulsion only 
(ix, vortex — 2 X 10^'*P_ii erg s^^) with the source con- 
tinuing to evolve along the vortex band (see Figure 2 in 
OLNIII). On the other hand, a typical AXP (older ring- 
bearing source) quiescent phase is dominated by emission 
from the HS induced by accretion from the outer edge of 
the ring with Lace, given by eq.®. These will evolve hor- 
izontally (at constant Lhs ^ W'^^ erg s^^) as discussed in 
OLNII (see Figure 2 in OLNIII). 

It is only during the bursting phase that transient and 
regular AXPs would look the same since they both expe- 
rience the BF effect and the related accretion and feed- 
back process. In the transient case, as the BF dies out 
when reaching the outer edge, the system becomes dom- 
inated again by vortex expulsion while regular AXPs re- 
sume their accretion dominated (from the outer edge) qui- 
escent phase. 

6.2. Birthrate 

In our model, a QS-ring system experiences an X- 
ray/radio outburst every few hundred years. Since we have 
observed 2 in a few years, located at distances of 3 to 9 kpc 
away (meaning we see about 30% of them), it is suggestive 
of a rate of 1 per year for the whole galaxy. This implies a 
total population of (1/yr) x 100 yrs ^ 100 in the galaxy. 
Since the ring will be consumed on a timescale of ^ W^- 
10^ years (see eq. 5 in OLNIII) this gives a birthrate of (1 
per 1000 years) to (1 per 100 years). Within uncertainties, 
the birth rate of transient AXPs derived above is consis- 
tent with the expected birthrate of AXPs (~ 1/300 yrs; 
Gill&Heyl 2007; Leahy&Ouyed 2009). As we argued in the 
previous section, in our model transient AXPs evolve into 
typical AXPs thus sharing the same birthrate. 

6.3. Model features and predictions 

The general predictions, starting with the X-ray emission, 
in our model are: 

- The 2 Blackbodies : 

^ Inside the star, the field lines expand outwards following 
the vortices. 
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Overall, during burst and quiescence, the cool BB 
(from Keplerian ring; or from the co-rotating shell as 
in OLNI) arises from reprocessing radiation of the hot 
BB from the star (either vortex annihilation or ac- 
cretion HS). We suggest, the relations in Figure 6 in 
Nakagawa et al. (2009) can be explained in the con- 
text of reprocessing during both quiescence and burst. 
Table 1, summarizes the different emission components 
for different objects in different states in our model. 
Compared to SGRs (see OLNI) and transient AXPs, 
typical AXPs acquire an additional BB from the HS 
during quiecence. Only during bursting do transient 
AXPs acquire a HS. As for SGRs in our model, we 
recall (see OLNI) they are born with a co-rotating 
shell (i.e. non-Keplerian degenerate ring; see OLNI). 
Future work will consider hard X-ray emission from 
non-thermal processes related to magnetic reconnec- 
tion following vortex expulsion and to accretion onto 
the QS. 

— During X-ray Burst: Since the accreted material in our 
model consists mostly of dissociated iron {Z ^ 13), 
we predict some sort of signatures either during chan- 
neling along the field line or on impact on the HS 
- maybe absorption lines or proton cyclotron lines 

0.14i?Qg,i4 keV) from any of element in the Ne- 
to-S group (e.g. Ne, Al, Si). These signatures should 
be common for both transient and typical AXPs in our 
model. 

— During X-ray Quiescence: When edge accretion occurs 
the signatures (e.g. absorption lines or proton cy- 
clotron lines) should be from the Z — 26 rather than 
the Z ~ 13 nuclei. 

— X-ray variability : During both the quiescent and 
bursting phases. X-rays from the ring-|-atmosphere sys- 
tem should be unpulsed and may carry the Keplerian 
sgnature via millisecond variability. 

The general predictions for the radio emission are: 

— Radio Delay: As can be seen from equation p^ . the 
shorter the period of the star the smaller the delay, 
idciay, between the X-ray burst and the radio emis- 
sion. Interestingly, the 2 radio emitting AXPs so far 
are observed are those with the smallest period. 

— Radio Flux : Furthermore, as can be seen from cq. p7|) . 
as P gets larger the radio following X-ray burst gets 
very faint making it more difficult to detect. Combined 
with long fdciay for large P, we argue these to be 
the reasons why radio is not observed following X-ray 
bursts in AXPs with higher (^ 10 s) period. 

Finally, we list specific predictions in our model in the 
case of XTE J1810-197: 

— XMM warm BB : The XMM warm BB for XTE 
J1810— 197 (ring+atmosphere system in our model) 
will evolve back to one single BB (ROSAT BB) with 
temperature Tq. 



Table 1. Thermal components in our model 



Sources 


Quiescent phase 


Bursting phase 


SGRs 


2 BBs 


3 BBs 


Transient AXPs 


2 BBs 


3 BBs 


Typical AXPs 


3 BBs 


3 BBs 



- XMM hard BB : The XMM hard BB (the accretion 
HS in our model) will disappear following burst (i.e. 
once the Bohm front reaches i?out)- However, we ex- 
pect edge effects to appear (e.g. flattening or even jump 
in the hot BB area; see last panel in Figure 3) just be- 
fore accretion shuts off. 



6.4. Further implications 

There exist two aspects of our model that might provide 
some answers to fundamental issues in pulsar magneto- 
spherefl First, the source in the quark-nova model is born 
as an aligned rotator and secondly is the fact that here the 
degenerate ring, by Keplerian shear, is a natural source of 
plasma (carried by the bubbles to the I.e.) for the magne- 
tosphere. Also, pairs are naturally supplied to the magne- 
tosphere by vortex annihilation. This is left as an avenue 
for future investigation. 

6.5. Model Limitations 

Our model suffers from a few caveats : (i) The interac- 
tion between the BF and the magnetic field in the atmo- 
sphere is rudimentary at this stage. Understanding the 
exact mechanism of feeding material from degenerate ring 
into the atmosphere and onto the field lines within the 
Bohm front is essential as this process is behind the feed- 
back process between the ring and the HS on the star. This 
is a complex problem/system that would require detailed 
MHD simulations before we can confirm this aspect of 
our model; (ii) While the bubble generation mechanism 
we propose is common in systems involving dipole files 
threaded by a disk on one end and a star on the other, 
the outward propagation of the bubbles within the I.e. de- 
pends crucially on conditions within the magnetosphere. 
In particular, the value a = 1 was chosen solely on the 
assumption that the bubbles will propagate in an equa- 
torial wind (given the Keplerian nature of the bubbles at 
birth) within the magnetosphere. For highly magnetized, 
aligned rotators, the physics of the magnetosphere within 

^ A disk model for X-ray emission from pulsars was consid- 
ered by Michel & Dessler (1981). However, their disk is funda- 
mentally different from our ring in that, where they hypothe- 
size electron degenerate material (left over from a supernova), 
in our model the material is completely relativistic-degenerate 
(from a Quark- Nova). This difference has many consequences, 
the most notable are a more efficient accretion mechanism 
(conversion energy from the hadron to quark transition as 
well as gravitational energy is released), and a slower (viscous) 
spreading rate for the ring. 
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the Lc. is yet to be understood and the problem solved. 
Furthermore, the bubble dissipation mechanism as they 
reach the I.e. remains to be demonstrated. The detailed 
structure of the thin shear layer at the light cylinder is 
not well studied and would require numerical simulations. 
For now our explanation is based on the assumption that 
the thin transition shear layer at the I.e. slowly destroys 
the bubble by shearing off the part that touched the light 
cylinder. 

7. Conclusion 

There are two fundamental components in our model 
for AXPs and transient AXPs namely, the QS and the 
Keplerian ring. In quiescence, vortex annihilation on the 
QS gives rise to thermal and non-thermal X-ray emission. 
The ring reprocesses the emission to give a second cooler 
BB emission component. Outburst is triggered by accre- 
tion of a small inner part of the ring (i.e. the wall). The 
two main consequences are production of light (Z ~ 13) 
nuclei and triggering MHD accretion onto the QS (yield- 
ing the HS). The interplay between the Bohm diffusion 
(i.e. i?BF term) and depletion of light nuclei (i.e. fi) gives 
rise to a rich behavior, necessary in order to account for 
the observed behavior of XTE J1810-197. Finally, one 
can ask if such a small Keplerian degenerate iron-rich ring 
could form around a neutron star. Ring formation when 
the neutron star is born appears implausible since a proto- 
neutron star is large compared to the ring size. After for- 
mation, there is no obvious mechanism to eject degenerate 
material unless a violent change of state, like a QN occurs. 
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Appendix A: The Bohm diffusion front 

The hot front propagates from the inner parts of the ring 
at a Bohm diffusion rate given by 



Ar = Vbf 6.25 x 10 



— 4 ^kcv 



1/2 



t 



1/2 



(A.l) 



where B±,i3 is the magnetic field component perpendicu- 
lar to the front motion in units of 10^^ G. The parameter 
ipBF carries some uncertainty related to the fact that the 
above formula is semiempirical (i.e. the constant 1/16 in 
the original Bohm diffusion expression has no theoretical 
justification; see §5.10 in Chen 1984). This parameter is 
derived self-consistently in our model as explained in 
below. For a qualitative derivation of Bohm diffusion we 
refer the interested reader to Alexeff&Rader (1991) and 
Belyaev (2001). 

The ring geometry allows us to write i?_L,i3 — 
tan (ac)-BBF,i3 where to a first approximation tan (etc) — 

(-ffftm.,out - -f^atm.,in)/(^out " -Rin) and BbF ^ -Bav. = 

5Qs(-RQs/^av.)^ where i?av. = (i?i„+i?out)/2; hereBqs = 



VSuPP and Rqs are the star's magnetic field and radius, 
respectively. The temperature cancels out from eg. dA.ll) 
since B± oc iJft^i Ty^cv/fJ- so that. 



1/2 



Ar 10"^ km tAbf 



A*3.3 



B: 



av.,13-Rout,15 / 



1/2 



1/2 

days 



(A.2) 



Appendix B: Centrifugal ejection and channeled 
accretion 

As the system (ring-|-atmosphere) is heated, the boundary 
between the non-degenerate atmosphere and the degen- 
erate Z 13 layer moves downwards into higher density 
layers (see eq.((T])). Thus Keplerian ring material is fed into 
the atmosphere and is ejected as a magnetohydrodynamic 
(MHD) wind as outlined below (see also lower panel in 
Figure 2). 

At keV temperatures the gas is sufficiently ionized 
everywhere in the atmosphere that ideal MHD can be 
used. Since the Lorentz force only has components per- 
pendicular to the field, the gas is free to move along the 
co-rotating field line under the influence of other forces. 
Under these conditions, it has been shown that the wind 
can be launched centrifugally if the field direction is in- 
clined at an angle less than 60° to the radial direction 
(Blandford&Payne 1982). Using conservation of specific 
angular momentum (oc r^fi), a non-degenerate Keplerian 
particle (i.e. fix) loaded at a footpoint ro will be flung 
out to larger radii; the so-called "bead-on-wire" analogy. 
These conditions are easily met in our model so that the 
newly unveiled Keplerian material (following heating of 
the degenerate ring) finds itself threaded by the highly in- 
clined magnetic field and is fiung out centrifugally. The 
wind is then channeled towards the star's surface by the 
strong dipole. 

The mass flux, rriwind, is regulated by conditions (i.e. 
density) at the slow mode critical point: 



47ri?BFCsp(Patm.,b-ffi 



atm. ,b 



-Patm.,q-ffatm.,q)^atm,th , (B.l) 



where (/3atm.,b-ffatm.,b — Patm.,q-ffatm.,q) = Patm. ,b^^atm. ,b X 

/K(7b) describes the mass injected from the Keplerian ring 
into the atmosphere. The function /K(7b) = l—Tq^^/T^^"^ 
evolves from ~ 0.8 to zero as the non-degenerate Keplerian 
material is depleted and the system evolves back to its 
quiescent phase. In a sense it acts as a "valve" that shuts- 
off accretion after the heating front has passed. 

rriwind, can be recast to 



The accretion rate, rhacc. 



1.77 X 10i^Csp,o.oi 



,5/2 
BF,15 



3/2 



(B.2) 



The factor carries the uncertainty on the exact loca- 
tion of the slow MHD mode. We will adopt ^sp ~ 0.01 
as inferred from 2-dimensional (Ouyed et al. 1997) and 
3-dimensional (Ouyed et al. 2003) MHD simulations of 
disk winds. The slow MHD point governs the mass loss 
rate. The degenerate ring responds by suppling mass at a 
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rate ttibf /3atm27ri?BF'^BF^^i where v± is the velocity of 
degenerate Keplerian material crossing the boundary to 
non-degeneracjH, and dep is the radial width of the BF 
as illustrated in Figure 2. 

B.l. Accretion shut-off 

The accretion relies upon new Keplerian material being 
fed from the degenerate ring to the atmosphere in the BF 
region, and ends once the magnetic field has re-inforced 
co-rotation of the atmosphere. This defines the inner edge 
of the BF as illustrated in the lower panel of figure 2. The 
penetration timescale can be estimated from equation (16) 

in OLNII to be tb,k ~ 4 days i?K.cm/'?keV where the depth 
of the non-degenerate Keplerian layer Hk ^ lO^^Hi^ ~ 
10~^7?ring is in units of centimeters. The size (i.e. radial 
width) of the BF front is found from dep — '^bp'''b,k ^ 10^ 
cm which is of the same order as the depth of the Z ~ 13 
layer. 



B.2. Depletion timescale of the Z ~ 13 nuclei 

We recall that following wall accretion and subsequent 
irradiation of the ring, up to 10'*^ of iron nuclei are dis- 
sociated forming the Z ^ 13 layer. Subsequent depletion 
of the Z ^ 13 nuclei by accretion leads to a return to an 
iron-rich atmosphere (i.e. /Zq ~ 3.3) as shown in OLNII: 



1 1.1 Is f/. 

- = — + ( )exp-*/^i^ 

M Mq Mb Mq 



(B.3) 



where T13 is the depletion timescale given by ri3 ~ 
wia/macc. where mis ~ N13 x (13 mn). For typical val- 
ues following wall accretion /i = /ib ~ 2.3, i?BP ^ Rim 
Ybf ^ 0.4, and ij — 0.1, we find 



ri3 - 1.4 X 10"* days 



13,46 



(B.4) 



■sp,0.01 



where -/Vi3,46 is the total number of dissociated iron nuclei 
in units of 10^^. 

We have already estimated the size of the BF region 
above which we found to be of the same order as the depth 
of the Z ~ 13 layer which implies that our scenario self- 
consistently leads to the accretion of most of the Z ~ 13 
material during burst (i.e. during the time it takes the BF 
to comb the ring's surface). In other words, the depletion 
timescale can also be estimated from eq. (|X2t by taking 
Ar = (i?,out — R'm) and /i = /^b ^ 2.3 to get 



Tl3 



1.5 X 10^ 

^Ip 



days (i?o 



R 



inJkm-fi^out,15^av.,13 j 



1/2 



(B.5) 



® We recall that in the degenerate Keplerian disk, lower vis- 
cosity and conservation of angular momentum imply mass flows 
in the ring is along lines of constant angular momentum which 
is nearly vertical (i.e. z-direction; see Appendix A in OLNII). 
Thus as the non-degenerate Keplerian material is "sucked up" 
in the wind, more material is supplied almost vertically from 
the underlying Keplerian material. 



where (i?out ~ ^in) is in units of kilometers. Equating the 
previous two equations we find 



C^p , 0. 1 ( ^out ^in ) km ^out , 1 5 ^av . ,13 



BP 



1/2 
13,46 



(B.6) 



which fine tunes the speed of the BF in our model (i.e. 
ea.(IX2])). 

Appendix C: Ring-Pole interaction: the hot spot 

As illustrated in Figure 2, there is a direct link between 
the BF moving outward along the ring surface (increasing 
its area) while decreasing the area of the HS on the surface 
of the star. It is straightforward to show that the total HS 
area (both poles) is 



A 



HS 



47ri?Qg X (cos/3bp 



cos An) , 



(C.l) 



where the angle /3 defines the colatitude (i.e. measured 
from the polar axis) on the star's surface where material 
channeled along a given field line lands . We approximate 
the field by a dipole configuration, giving us 



1/2 



smof 



where 



^ , -^ring 

a = arctan — — — 

2 R 



(C.2) 



(C.3) 



is the angle from the polar axis to the footpoint, at radial 
distance R, of the field line where the Keplerian material 
is loaded; -ffring ^ 2.68 kmp^^f^ 9^15^ is the ring's vertical 
height as given in H2.'2[ For our fiducial values, i?in = 15 
km, i?out = 25 km and Rqs = 10 km, we find (cos /3out — 
cos/3in) '~ 0.192 which gives an initial area of the HS at 
one pole of the order of ^ 120 km^. 
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